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Abstract 

A study is carried out of the role of the ahgned neutron-proton pair with angular momentum 
J = 9 and isospin T = in the low-energy spectroscopy of the N = Z nuclei ^^Cd, ^^Ag, and 
^^Pd. Shell- model wave functions resulting from realistic interactions are analyzed in terms of a 
variety of two-nucleon pairs corresponding to different choices of their coupled angular momentum 
J and isospin T. The analysis is performed exactly for four holes (^^Cd) and carried further for 
six and eight holes (^^Ag and ^^Pd) by means of a mapping to an appropriate version of the 
interacting boson model. The study allows the identification of the strengths and deficiencies of 
the aligned-pair approximation. 

PACS numbers: 21.60.Cs, 21.60.Ev 
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I. INTRODUCTION 



The study of nuclei with equal numbers of neutrons and protons (A^ = Z) is one of the 
declared goals of radioactive-ion-beam facilities, either in operation or under construction. 
Such nuclei are well known when and Z are small but, as the atomic mass number 
A = N+Z increases, they lie increasingly closer to the proton drip line and therefore become 
more difficult to study experimentally. Nevertheless, several phenomena of interest, such as 
the breaking of isospin symmetry or the emergence of new collective modes of excitation, 
are predicted to become more pronounced with increasing A, and this constitutes the main 
argument for undertaking the difficult studies of ever more heavier N = Z nuclei. 

Arguably, the goal of most interest in this quest is the uncovering of effects due to isoscalar 
(T = 0) neutron-proton (n-p) pairing. In contrast to the usual isovector (T = 1) pairing, 
where the orbital angular momenta and the spins of two nucleons are both antiparallel {i.e., 
L = and 5 = 0), isoscalar pairing requires the spins of the nucleons to be parallel {S = 1), 
resulting in a total angular momentum J = 1. Collective correlation effects are predicted to 
occur as a result of isoscalar n-p pairing but have resisted so far experimental confirmation 
because (i) the jj coupling scheme which is applicable in all but the lightest nuclei disfavors 
the formation of isoscalar n-p pairs with L = {2!, and (ii) the states associated with this 
collective mode are of low angular momentum J and often hidden among high-J isomeric 
states which hinders their experimental detection. 

Currently, N = Z experiments are approaching ^''^Sn, involving studies of nuclei such 



as Pd [3[ where nucleons are dominantly confined to the lgg/2 orbit. In the context of 
these experiments, Blomqvist recently proposed ^ that a realistic description of shell-model 
wave functions can be obtained in terms of isoscalar n-p pairs which are completely aligned 
in angular momentum, that is, pairs with J = 9. This proposal is attractive since it is 
well adapted to the jj coupling scheme, valid in this mass region, and because it should 
encompass the description of high-J isomeric states. Blomqvist's idea is related to the so- 
called stretch scheme which was advocated a long time ago by Danos and Gillet jsl. In 
the stretch scheme, shell-model states are constructed from aligned n-p pairs, treated in a 
quasi-boson approximation which neglects antisymmetry between the nucleons in different 
pairs. The latter approximation is absent from Blomqvist's approach. 

In this paper we examine Blomqvist's proposal for the N = Z nuclei ^^Cd, ^^Ag, and 
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^^Pd. We consider several realistic two-body interactions for the lgg/2 orbit and analyze 
the shell-model wave functions, obtained with these interactions, of the four-nucleon-hole 
system (^^Cd), in terms of a variety two-pair states. For the six- and eight-hole nuclei (^^Ag 
and ^^Pd) a direct shell-model analysis in terms of pair states is more difficult, and we prefer 
therefore to carry out an indirect check by means of a mapping to a corresponding boson 
model. In these cases our approach is intermediate between that of Blomqvist |^ and of 
Danos and Gillet jsl. The boson mapping takes care of antisymmetry effects in an exact 
manner on the level of four nucleons but becomes approximate for more. 

This paper is organized as follows. First, some necessary concepts and techniques are 
introduced: the formulas needed to carry out a shell-model calculation in a pair basis are 
given in Sect. Hi] and two mapping techniques from an interacting fermion to an interacting 
boson model are reviewed in Sect. IIIII The results of our analysis of N = Z nuclei are 
presented and discussed in Sect HVl The conclusions and outlook of this work are summarized 
in Sect. |Vl 



II. FOUR-PARTICLE MATRIX ELEMENTS 



This section summarizes the necessary ingredients to carry out calculations in an isospin 
formalism. The formulas given are valid for fermions as well as for bosons. Four-particle 
states are described by grouping the particles in two pairs. These two-pair states can be 
used, for example, as a basis in a shell-model calculation, facilitating the subsequent analysis 
of the pair structure of the eigenstates. Furthermore, the two-pair representation of four- 
particle states is the natural basis to map the shell model onto a corresponding model in 
terms of bosons. Once this mapping is carried out, the original interacting fermion problem 
is reduced to one of interacting bosons, which can also be solved with the results summarized 
in this section. 

In the pair representation of four particles with angular momentum j and isospin t 
(both integer for bosons and half-odd- integer for fermions), a state can be written as 
|(jt)^(JiTi)(jt)^(J2T2); JT) where particles 1 and 2 are coupled to angular momentum and 
isospin JiTi, particles 3 and 4 to J2T2, and the intermediate quantum numbers JiTi and 
J2T2 to total JT. It is convenient to introduce a short-hand notation that takes care 
simultaneously of angular momentum and isospin quantum numbers and we denote JT 
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as T; indices will always be carried over consistently, i.e., Fj refers to JjTj, 7 to jt, etc. 
The state \{jtf{JiTi){jtY{J2T2)]JT) is then denoted as \-f'^{Ti)-^'^{T2)-J) . This state is 
not (anti) symmetric in all four bosons (fermions) but can be made so by applying the 
(anti)symmetrization operator P, 

= E[7'(ra)7'(r.);r|}7'[rir2]r]|7^(r„)7^(r,);r), (i) 

where Y^'^(Va)l'^{Xh)', r|}7''[rir2]r] is a four-to-two coefficient of fractional parentage (CFP). 
The notation in square brackets [rir2] implies that ([1]) is constructed from a parent state |6| 
with intermediate angular momenta and isospins JiTi and 72^2- It is implicitly assumed 
that all pair quantum numbers are allowed, i.e. that Jj is even for Tj = 1 and odd for Tj = 0. 
The four-to-two CFP is given by 



[7^(r,)7^(r,);r|}7lrir2]r] 
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where a = +1 for bosons and a = —1 for fermions. Furthermore, A/'^p^^ is a normalization 
constant, ^r^Ta = ^JiJa^TiTa, ^Vir, = ^JiJ^TiU, and 
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where the symbol in curly brackets is a nine-j symbol and [Fj 
normalization constant is known in closed form as 



(2Ji + 1)(2T, + 1). The 
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l + (-)''+^5r,r,+4a 
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The states ([1]) do not form an orthonormal basis. One needs therefore to determine the 
overlap matrix elements which can be written in terms of the CFPs as 



(7^[FiF2]F|7^[F;F'2]F) = Y: [7'(ra)7'(r.);F|}7lFiF2]F][7^(F,)7^(F,);F|}7lF;F^]F] 



Furthermore, the matrix elements of the two-body part of the Hamiltonian can be expressed 

as 

TaA 

(6) 

where z/a are two-body matrix elements between normalized two-particle states, z^a = 
(7^; A|if2|7^; A). The label A is a short-hand notation for the two-particles' angular mo- 
mentum A and isospin T. For example, in the nuclear shell model the particles are nucleons 
with half-odd- integer angular momentum j and isospin t = 1/2. We recall that in this 
case the antisymmetric two-nucleon states are uniquely determined by the total angular 
momentum A and that the total isospin is a redundant quantum number which is T = for 
odd A and T = 1 for even A. The two-body matrix elements therefore depend on A only, 
ux = (j^; X\H2\j'^', A). In the interacting boson model, in contrast, the particles are bosons 
with integer angular momentum j and integer isospin t = or t = 1. In the latter case, 
the two-particles' isospin T (obtained by coupling t = 1 with t = 1) is not redundant but is 
needed to fully characterize the two-particle state. 



III. METHODS OF MAPPING 



The basic and common idea of different boson mapping methods is to truncate the full 
shell-model space to a subspace which is written in terms of fermion pairs and to estab- 
lish subsequently a correspondence between the fermion-pair space and an analogous space 



which is written in terms of bosons. This section reca^ 
OAI (Otsuka-Arima-Iachello) [3] and the democratic 
applications. 



Is briefly two methods, namely the 
] mappings that will be used in the 



A. The OAI mapping 

We limit ourselves here to the problem of n nucleons with isospin t = 1 /2 in a single orbit 
j which is the case of interest in this paper. A detailed description of the method in its full 
generality can be found in Ref. 



Let us first introduce the pair creation operators 

AiJT) _ t UiJT) /^x 

^MjMt ~ ^2 \MjMt^ v ) 

where is the creation operator of a t = 1/2 nucleon in orbit j and the square brackets 
denote couphng to a tensor with angular momentum J and isospin T, and projections Mj 
and Mti respectively. For nucleons in a single orbit j, the pair (JTj) is totally determined by 
its angular momentum J and its isospin follows from it (z.e., T = for odd J and T = 1 
for even J). The starting point of the OAI mapping is a given shell-model Hamiltonian 
which is a scalar in J and T, and the selection of a number of fermion pairs Ji, J2, ■ ■ ■ , Jp, 
to each of which will be associated a boson with the same angular momentum Jj. (For 
the bosons this angular momentum is by convention denoted as £j, £i = Jj.) For a single-j 
orbit one may ignore the single-particle term of the Hamiltonian which we therefore take 
to be of pure two-body character and denote it as i^f . The one-body term of the mapped 
boson Hamiltonian is directly obtained from the matrix elements of the shell-model 
Hamiltonian H2 between the pair states, 

e,^ = {i^\Hy.) = {Ji\H!\Ji). (8) 

These matrix elements are interpreted as single-boson energies. 

The determination of higher-body terms of the mapped boson Hamiltonian is more com- 
plicated, as we illustrate here for the two-body part. For a given angular momentum J 
and isospin T, one first enumerates all possible two-pair states \Fi) = \ JaJb] JT), where the 
index i = 1, . . . ,d corresponds to pairs (a, b) with 1 < a < b < p. These will be mapped 
onto the two-boson states \Bi) = \iah'i JT), 

\Fi)^\B,), t = l,...,d. (9) 

This correspondence cannot be established directly since the boson states are orthogonal 
while the fermion states are not: 

{Bi\Bj) = 5^j, Oij = {Fi\F,) ^ 5ij, (10) 

where Oij is the overlap matrix in fermion space. To deal with this problem, the strategy 
of the OAI mapping is to define for each JT a hierarchy of pair states, to apply a Gram- 
Schmidt orthogonalization of this ordered sequence, and to associate the boson states with 
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the orthonormalized fermion states. This leads to the definition of the correspondence 



i = l,...,d, (11) 
where the states form an orthonormahzed basis. One has the series 



\F2) = , J \F2)-d2,\F,)), 

yo22 - 



\Fk) = ^^[\Fk)-'^^dk.\F,)], (12) 



until k — p, with 



fc-i 



= Okk - ^ki = {Fk\Fi). (13) 



1=1 



An efficient algorithm to carry out the Gram-Schmidt procedure involves expanding the 
orthonormal states in the non-orthogonal basis as 

k 

\Fk)^J2<'ki\Fi), (14) 

i=l 

and calculating the coefficients a^i (needed for k >i) recursively from 

2 -j^ fc— 1 i 

0-kk — rr-ri flfei — 7x7= X] ^ki'fJ'i'i, Ofej = ^ On'Ofej', (15) 

which as only input requires the knowledge of the overlap matrix elements Oij. The matrix 
elements of the fermion Hamiltonian in the orthonormal basis \Fk) are then given by 



{Fk\H!\Fi) = J2(ikiaij{Fi\Hi\Fj), (16) 

ij 

and the mapped boson Hamiltonian follows from 

(44; JT|^^|4'4'; JT) = {JJb] JT\Hl\jyjv- JT). (17) 

Note that this equation defines the entire boson Hamiltonian up to and including two-body 
interactions. To isolate its two-body part one should subtract the previously determined 
one-body terms according to 

(44; JT\H^\L'h'\ JT) = (44; JT\H^\L'h'-, JT) - (e,„ + e,;)5^,5w, (18) 
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always assuming that a <b and a' <b'. 

This is the version of the OAI mapping as it will be applied in this paper. The boson 
Hamiltonian is obtained from the two- and four-particle systems and is kept constant for 
systems with higher numbers of particles. This is akin to the democratic mapping but the 
latter has the advantage that no hierarchy of two-pair states is required, as is discussed in 
the next subsection. 

Nevertheless, the hierarchy imposed in the OAI mapping is often based on arguments 
of seniority which enables an extension to n-particle systems. To illustrate this point, we 
consider the case of a truncation onto a subspace constructed out of pairs of angular momenta 
J = {S pair) and J = 2 {D pair). (In this example we assume for simplicity identical 
nucleons so that isospin can be omitted.) The creation operators 5*^ and D'^ are defined as 

S^ = aT\ D^m = A'S'. (19) 
The truncated shell-model space is the SD subspace spanned by the states 

IjX; JMj) = -±p(5t)("-'')/2[(Z}t)-/2](2|o), (20) 

where |o) refers to the closed core, ^ denotes additional quantum numbers related to the 
intermediate angular-momentum couplings of the D pairs, and A/p is a normalization con- 
stant. The operator V is needed to ensure orthogonality of the basis. Its effect is a projection 
onto a subspace with seniority v, the number of particles not in pairs coupled to J = 0. 
As a result the states fl20l) are normalized and can be mapped onto sd states. This can be 
expressed by the general correspondence 

IjX; JMj) ^ In.n.C; JMj) = -l(st)"=[(cit)«<^]g|o), (21) 

where |o) is the boson vacuum, Ug and are the s- and cZ-boson numbers, respectively, with 
Us = {n — v)/2 and rid = and A/r is a normalization constant. The matrix elements of 
the mapped boson Hamiltonian are now defined as 

{n^ridi; JMj\H''\nWd'- JMj) = (jX; JMjlH^lfv'^'; JMj). (22) 

The difference with Eq. ( ITTl) is that the latter equation is specific to n = A particles while 
the mapping ( 122|) applies to any n. The calculation of an n-particle fermion matrix element 
is possible because the seniority formalism allows it to be reduced to n = max(w,w') 
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Hence, as long as v,v' < 4 the fermion matrix element can be computed and fixes the 
corresponding boson matrix element. The advantage of this procedure is that it determines 
a boson Hamiltonian which varies with particle number yielding a more adequate description 
of Pauli effects. Its disadvantage is that a generalization towards any choice of pairs and/or 
to neutrons and protons is difficult, if not impossible. Therefore, we opt in this paper for 
the simpler OAI mapping that defines a constant boson Hamiltonian from the two- and 
four-particle systems. 



B. The democratic mapping 

Again the starting point is the selection of a number of fermion pairs Ji, J2, ■ ■ ■ , Jp and a 
corresponding series of bosons with energies eg. determined from Eq. ([HD. To establish the 
correspondence complicated by the non-orthogonality of the fermion two-pair states, 
the democratic mapping relies on the diagonalization of the overlap matrix Oij defined in 
Eq. ffTOj) . This diagonalization provides, besides the eigenvalues Ok, an orthogonal basis \Xk), 

= ^Cfc,|F,), A; = l,...,d, (23) 

i 

with 

{Xk\Xi) = J2 CkiCijOij = OkSki, k,l = l,...,d. (24) 

We follow the convention of labeling non-orthogonal basis states by i or j, and orthogonal 
basis states by k or /. The coefficients Cki = {Fi\Xk) = {Xk\Fi) are transformation coefficients 
from the non-orthonormal basis \Fi) to the basis \Xk) and they satisfy 

^CkiCii = 5ki, k,l = l,...,d. (25) 

i 

The states \Xk) form an orthogonal but non-normalized basis since they satisfy (IMj) . To 
normalize these states, we define 

\Fk) = ^\Xk), k = l,...,d, (26) 



. Ok 

which indeed satisfy 

{Fk\Fi) = 5ki, fc,/ = l,...,rf. (27) 

The basis \Fk) has the same role as the basis \Fk) in the OAI mapping but it is important 
to realize that both bases are not identical and lead to different boson Hamiltonians. The 
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democratic mapping relies on the definition of a transformation in boson space which is 
analogous to the one in fermion space 

= 5]cfc,|5,), k = l,...,d. (28) 

i 

From the orthogonality of the basis \Bi) and the properties of the coefficients Cki it can be 
shown that these states form an orthonormal set, 

{Bk\Bi)=5ki, kj = l,...,d. (29) 

We have now arrived at fermion and boson bases that are both orthonormal, and we can 
therefore establish the mapping 

\Fk)^\Bk), k = l,...,d, (30) 

and determine the matrix elements of the boson Hamiltonian in this basis, 

{Bk\H''\Bi) = {Fk\H!\Fi), k,l = l,...,d. (31) 

With use of the inverse of the relation ( 125]) . of the equahty (IHTj) . and of Eqs. ( 12^ and 

the matrix elements of the boson Hamiltonian in the original basis can be written in terms 

of those of the fermion Hamiltonian (also in the original basis) as 

{Bi\H^\Bj) = / CkiCki'CijCij' HJ,-i , (32) 

kl i'f V^kOl 

where H^j, = {Fii\H2\Fji) . Again, as with the OAI mapping, this defines the entire boson 
Hamiltonian from which the two-body part can be isolated by applying Eq. (fT8|) . 

In some cases, m fermion vectors \Fi) are linearly dependent on the d — m others, leading 
to m vanishing eigenvalues of the overlap matrix, 

Oj = 0, i = 1, . . . ,m. (33) 

This problem can be solved by excluding from the fermion space m states and by 
calculating the matrix elements of from the remaining d — m states. All other matrix 
elements of are defined such that 

{Bi\ \ Bi) = oo, i = d-m + l,...,d, 

mH''\B,) = {B,\H''\Bi) = 0, J, (34) 

where i and j in the last equation take the values i = 1, . . . , and j = d — m + 1, . . . ,d. One 
has still to decide which m vectors to remove from the fermion space. To avoid arbitrariness, 
one chooses the m fermion states with the smallest overlap with the shell-model states. 
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IV. APPLICATION TO N = Z NUCLEI 



Of particular interest in this work are N = Z nuclei with neutrons and protons in the 
1(79/2 orbit. Recently, Blomqvist [4| has conjectured that a valid interpretation of yrast 
states in these nuclei can be obtained in terms of n-p T = pairs which are coupled to 
maximum angular momentum J = 9 and which therefore can be termed aligned n-p pairs. 
This is contrary to the usual interpretation of such states which involves low-spin pairs with 
isospin T = 1 and possibly also with T = 0. In this section we examine Blomqvist's proposal 
with specific reference to the nuclei ^^Cd, ^^Ag, and ^^Pd, corresponding to four, six, and 
eight holes with respect to the ^°°Sn core, respectively. Our study consists of two separate 
parts: the analysis of shell-model wave functions of ^®Cd in terms of a variety of two-fermion 
pairs and the mapping of shell-model onto corresponding boson states for ^®Cd and ^^Pd, 
following the formalism developed in the previous sections. First, an appropriate shell-model 
interaction should be determined. 



A. Shell-model interaction 

For the purpose of checking the stability of our results and the reliability of our con- 
clusions, we have carried out the analysis for three different shell-model interactions. The 
SLGTO interaction is taken from Serduke et al. 9l. It was used in the more recent analysis 

n 

of Herndl and Brown [10|] where it was found to give satisfactory results for the neutron- 
deficient nuclei in the mass region A = 86 to 100, which are of interest in the present 
study. A second shell-model interaction, which shall be named GF, is taken from Gross 



and Frenkel [llj. Both the SLGTO and GF interactions are defined in the 2pi/2 + 1(79/2 
shell-model space and, to carry out an analysis in terms of aligned pairs, it is necessary to 
renormalize them to the 1(79/2 orbit. The resulting two-body matrix elements are shown in 
columns 2 and 3 of Table [B The renormalization only affects the (A,T) = (0, 1) and (1,0) 
matrix elements since these are the only ones that also occur in the 2pi/2 orbit. 

To avoid the renormalization procedure, we define a third interaction directly for the 
1(79/2 orbit. This is done in the following way. The spectrum of ^°Nb is well known |l2| and 
enables the determination of the particle-hole interaction matrix elements 

uf = (1(79/2(1^79/2)-^ AT|4^|l(79/2(1^79/2)-^ AT). (35) 
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TABLE I 


Two-body matrix element 


s in the 


1(79/2 orbit, in units of MeV, derived 


irom the 


interactions SLGTO and GF, and from 


the experimental spectrum of ^'^Nb. 










= (ilanin)'^- XT\H^\{lQa,o)'^- XT) 




(A T) 


SLGTO 
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Nb90 




-2.392 




-2.321 


-1.758 




— 1.546 




— 1.524 


-1.225 




-0.906 




-0.937 


-0.573 




—0.747 




—0.700 


-0.521 


(4,1) 


—0.106 




—0.160 


0.064 


(5,0) 


-0.423 




-0.447 


-0.332 


(6,1) 


0.190 




0.140 


0.266 


(7,0) 


-0.648 




-0.640 


-0.481 


(8,1) 


0.321 




0.241 


0.334 


(9,0) 


-1.504 




-1.752 


-1.376 



The isospin T is determined from A, that is, T = 4 for all states except for A = which 
has T = 5. The absolute value of the matrix element fl35|) for the ground state (A = 8) is 
obtained from the binding energies of the surrounding nuclei as 



ph 



^(^"Zr) + E(^"Nb) - Ei^^m) - ErZi) 
-783.794 - 776.895 + 788.942 + 771.825 
0.078 MeV, 



(36) 



where E stands for the nuclear binding energy, taken from the 2003 atomic mass evalua- 



tion 



13| and corrected for the electrons' binding energy IJ]. (The minus sign is needed to 



convert from binding energy to interaction energy.) All levels A = 0, . . . , 9 of the particle- 
hole multiplet fl35|) are known in ^°Nb and fix the differences uf^ — u^^. A problem with 
this procedure concerns the choice of the relevant J'" = 1^ state in ^°Nb. Several of them 
are observed at low excitation energy and only one should be taken as a member of the 
particle-hole multiplet. We have taken the J'" = level at 2.126 MeV because this state is 



strongly populated in the /3-decay of ^°Mo (log ft = 4.9) [l5|. Furthermore, it appears to be 
the only 1+ state observed in the triton spectrum obtained in the charge-exchange reaction 
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90Zr(3He,t)90Nb Q. 

Once the particle-hole matrix elements are derived in this way, the particle-particle 
(or hole-hole) matrix elements are obtained from an inverted Pandya tranformation jl7 |. 
Since absolute (as opposed to relative) matrix elements have been extracted from binding en- 
ergies, care should be taken to use the appropriate Pandya transformation. Equation (18.63) 
of Ref. 6| gives the following relation between absolute matrix elements: 

j 3 A 



ph 



i?o-E(2A' + l){' ' y^r, (37) 
A' J J A 



where Eq is a constant {i.e., A- independent) interaction energy given by 

Eo = ^ E (2A' + iKf + 1^ E (2A' + iKf. (38) 

all A' J even A' 

To express the particle-particle in terms of the particle-hole matrix elements, as is needed 
here, the relation (137|) can now be inverted in the usual manner, leading to 

^A^^ = ^o-E(2A' + l){^ ^ ^K- (39) 
A' j j A' 

The problem with this relation is that it is of no use as long as the constant Eq is expressed 
in terms of particle-particle matrix elements, as in Eq. (|38ll . We need an equation for Eq in 
terms of particle-hole matrix elements. This can be obtained by inserting the expression f l5I?]) 
for in Eq. (155]) and solving for Eq, leading to 

Eq = y (2A' + l)z/P? - -^^^^i^t- (40) 

2j(2j + l)it' 2j(2j + l) ° ^ ^ 

The resulting particle-particle matrix elements are shown in column 4 of Table |T] under 
'NbQO', where the 'pp' index is omitted, ux = z/^^, as will be done from now on. Note that 
the matrix elements thus obtained differ by a constant from those of Sorlin and Porquet (see 
Fig. 6 of Ref. 18|), since we have taken care here of the constant interaction energy Eq. 

It is seen from Table [T] that the matrix elements of the SLGTO and GF interactions are 
similar. The biggest difference concerns the (A,T) = (9,0) matrix element which is more 
attractive by ~ 250 keV for GF. This conceivably might influence the approximation in 
terms of aligned A = 9 n-p pairs. The interaction derived from ^''Nb is less attractive (or 
more repulsive), in particular the pairing matrix element with (A,T) = (0,1). All three 
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TABLE II: Fractions of 2p]^/2 + ^99/2 shell-model eigenstates of the GF interaction that lie within 
the 1^9/2 subspace, expressed in percentages. 
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interactions, while having reasonable characteristics, are sufficiently different to test the 
robustness of our analysis. 

Before proceeding with the wave-function analysis, we first check to what extent the 
N = Z nuclei in this region can be described by confining nucleons to the lgg/2 orbit. This 
approximation should be reasonable for nuclei close to (south-west of) ^"^^Sn but will become 
increasingly poor as one approaches ^°Zr. The latter nucleus is known to be deformed 13] 
and hence a single (spherical) orbit will not suffice for a reliable description of nuclei in 
its neighborhood. To quantify the limitation to the lg9/2 orbit, we show in Table HT] the 
results of a shell-model calculation with the GF interaction. The table shows the fractions 
of 2pi/2 + 1(79/2 shell-model eigenstates of ^^Cd and ^^Pd that lie within the 1(^9/2 subspace, 
expressed in percentages. It is seen that for the lowest eigenstates this fraction is large, hence 
justifying the restriction to the 1(79/2 orbit. Therefore, we henceforth restrict the shell-model 
space to 1(79/2, in which case the formalism of Sect. [Ill can be applied. 

B. Shell-model analysis of the (159/2)^ system 

With the formalism developed in Sect. [Ill it is now possible to perform a shell- model 
calculation for four holes in the lgg/2 orbit and analyze the resulting wave functions in 
terms of two-pair states. We concentrate on the N = Z case which corresponds to two 
neutrons and two protons. For a given angular momentum J, the overlap matrix between 
all possible two-pair states |(jt)'^[(JiTi)(J2T2)] JT) is constructed using Eq. ([5]). The number 
d of linearly independent states JT is given by the number of non-zero eigenvalues of the 
overlap matrix. We then select d linearly independent but otherwise arbitrary two-pair 
states that shall be denoted in short as \Fi),i = 1,. . . ,d. Diagonalization of the overlap 
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angular momentum / 

FIG. 1: Overlaps of the (159/2)^ yrast eigenstates of the interactions SLGTO, GF, and Nb90 with 
angular momentum J and isospin T = 0, with the two-pair state J). 

matrix Ojj = (-Fj|Fj) allows the definition of an orthonormal basis 



(41) 



with Cfcj = Cfcj/ ^/ofc in the notation of Sect. Illli This basis can now be used to construct the 
energy matrix, 

kh\iil\Fi) = Y.-CkiCi^l,F,\El\F>), (42) 

to be computed with use of Eq. (E]) . The diagonalization of this matrix leads to the energy 
eigenvalues e^, r = 1, . . . , rf, with corresponding eigenvectors given by 



k ki 



(43) 



Note that the energies are independent of the choice of two-pair states \Fi) as long as 
the latter span the entire JT space. Although the above procedure might seem rather 
cumbersome for performing a four-hole shell-model calculation, it has the advantage of 
providing us directly with the pair structure of the shell-model eigenstates from the overlaps 



(^F{ I ) ^ ^ (^rk^kj Oij ■ 
kj 



(44) 



We now analyze the yrast eigenstates of the various interactions defined in Table HI that 
is, the quantities (FilEr)"^ for r = 1 and a variety of two-pair states Fj. Pairs with angular 
momentum J are generically denoted as Pj and explicitly, following standard spectroscopic 
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TABLE III: Coefficients flp;^jr2]r ™ expansion (I45p of the diagonal energies of the (9/2)'^ system 
for isospin T = 0. 
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notation, as 5, G, /, and K ioi J = 0, 2, 4, 6, and 8. Since this paper deals in particular 
with the aligned n-p pair with J = 9, we reserve for it the non-standard notation 'S' (for 
Blomqvist). The central result is shown in Fig. [1] which displays the quantity (Ji|i?^; J)^ 
where |Ji) is the yrast eigenstate with angular momentum J and isospin T = of the 
interactions SLGTO, GF, and Nb90. Most yrast states have a large overlap with i?^, as was 
shown by Blomqvist, but this is conspicuously not the case for J 8. It seems as if the 
two aligned n-p pairs do not like to couple to a total angular momentum which equals their 
individual spins. 

To acquire some insight in this finding, we note that simple expressions in terms of the 
two-body matrix elements are available for the diagonal energies of two-pair states, 

E{^'[T,T,]T) ^ {^'[T,T,]T\HimT,V,]T) = E4.r.]r ^a, (45) 

A 

where the coefficients oj^^p^jp of relevance to the present discussion are shown in Table lllli 
According to a recent paper by Talmi they are non-negative rational numbers; since 
they involve ratios of rather large integers in this case. Table IIIII gives numerical approxi- 
mations. Of all the coefficients given in this table, the important ones have A = 0, 1, and 
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TABLE IV: Overlaps of the (199/2)^ yrast eigenstates of the SLGTO interaction with angular 
momentum J and isospin T = with various two-pair states, expressed in percentages. 
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9, because these are the multipolarities of the most attractive interaction matrix elements. 
It is seen that the contributions of the aligned matrix element (A = 9) to the energies 
of the states | (9/2)^ [5^] J) and | (9/2) ^[SPj] J) remain more or less constant, independent 
of J . This is not the case for the pairing matrix element (A = 0) whose contribution 
to the energy of | (9/2)*^ [5^] J) disappears as J increases while it remains important for 
\{9/2Y[SPj\J) . The combined effect of these contributions is that the state | (9/2)^ [SPj] J) 
dips below |(9/2)^[52] j) around J ^ 8 and clS Qb result picks up the largest component of 
the yrast eigenstate. Hence, the feature of the disappearing B dominance around J ~ 8 is 
explained by a combination of geometry — the CFPs in the j = 9/2 orbit, and dynamics — the 
dependence of the interaction matrix elements on A. 

Figure [1] shows that the overlaps (Ji|i?^; J)^ are very similar for the three interactions. 
This finding is at the basis of the fact that the subsequent analysis gives consistent results 
for the three interactions. While there can be significant differences in the shell-model re- 
sults with the different interactions, the approximation in terms of aligned pairs is similar 
for the three interactions. In other words, if a particular shell-model state is well approxi- 
mated in terms of aligned pairs for one interaction, it is so for the other two as well; if the 
approximation is less good, it is so for all three interactions. Although we have carried out 
the complete analysis for the three interactions SLGTO, GF, and Nb90, we will show in the 
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following only the results of the former since it has a proven track record of satisfactorily 



io|. 



reproducing the data in the mass region of interest 

In Table HV] are shown the amplitudes in percentages for yrast eigenstates of the SLGTO 
interaction with even J and T = 0, that is, the quantities 100 x [Fi\Ej.Y for r = 1 and a 
variety of pair states The numbers illustrate that, at least at low and at high angular 
momentum J, the overlaps of the physical eigenstates with \B'^\J) are more important than 
those with other pair combinations. The percentages shown in Table IIVI also illustrate the 
non-orthogonality of the two-pair basis. For the example, the J = ground state has a 91% 
overlap with B"^ but also a 80% overlap with 5*^; this can only be if the overlap {B'^\S'^) 
itself is rather large. 



C. Boson mappings 

Ideally, one would like to perform a similar analysis of shell-model eigenstates for more 

enging problem, however, since it requires the formula- 



than four nucleons. That is a chal 
tion of a nucleon-pair shell model 
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22l | in an isospin-invariant formalism. In this paper 
we choose to extend our analysis toward higher hole number through the boson mapping 
techniques explained in Sect. Illli It is important to stress that this approximation goes 
beyond the original proposal of Blomqvist since it involves an additional assumption of the 
boson character of the fermion pairs. The results presented in this subsection therefore do 
not directly address Blomqvist's conjecture. 

Once the mapping is carried out for the two- and four-hole systems according to one of 
the two procedures described in Sect. IIIIl a Hamiltonian is obtained in terms of the selected 
bosons ii,i2, ■ ■ ■ Api which can then be applied to systems with two or more bosons. Such 
a description shall be referred to as C.1C.2 ■ ■ ■ ^p-IBM, where IBM stands for interacting boson 



model 23j. Note that all versions of IBM thus obtained are isospin invariant; for example, 
the sd-lBM is in fact the IBM-3 of Elliott and White Q. 

To compare the merits of different selections of fermion pairs, the following mappings are 
considered: 

1. A single fermion pair B with J = 9, T = 0, leading to the 6-IBM. 

2. Two fermion pairs S and B with J = 0,T = 1 and J = 9,T = 0, leading to the 



18 



TABLE V: Energies (in MeV) of T = levels for four nucleons in the lgg/2 orbit (^^Cd) calculated 
with the shell-model interaction SLGTO and compared with various versions of IBM obtained by 



democratic 
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ground state. 
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4.730 
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sb-lBM. 

3. Two fermion pairs S and D with J = and 2, both with T = 1, leading to the sd-IBM. 

4. Three fermion pairs 5", D, and G with J = 0, 2, and 4, all with T = 1, leading to the 
sdg-lBM. 

The first two cases are inspired by Blomqvist's conjecture, involving the aligned n-p pair, 
while the next two are the standard choice of the IBM |23i] and its most frequently used 



extension which includes g bosons. (For a review on the latter, see Ref. [24|.) 

The results obtained with the various boson Hamiltonians are compared with T = 
eigenstates of the SLGTO interaction for four, six, and eight nucleons in Tables |Vl IVll 
and IVIIl respectively. The numerical calculations have been performed with the codes 
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ArbModel [26| and IBM-3 [27]. The former is a general purpose program that can handle 
systems of fermions and/or bosons with arbitrary spins and can thus be used for the shell- 
model as well as the IBM calculations; the latter code is specifically written for the isospin- 
invariant sd-IBM. Alternatively, for three and four identical bosons {i.e., for 6-IBM) the 
calculations can be performed with the expressions given in Sect. [ITland equivalent ones for 
the three-hole case. 

A few remarks are in order. All results concern absolute energies. In the first line of 
each table are given the binding energies Eo of the T = ground states, as obtained in the 
various mappings, which should be compared with the corresponding quantity in the shell 
model. In subsequent lines are given the energies of a selected number of states, relative to 
this ground state. This might lead to some seemingly counterintuitive results. For example, 
it is seen from Table |IV] that the four-hole 2f shell-model state overlaps 97 % with a B'^ 
configuration. Why, then, should its excitation energy come out rather poorly in 6-IBM 
(0.678 MeV compared with 0.963 MeV in the shell model, see Table |V|) ? The reason is that 
the absolute energy of the 2i state is rather well reproduced (it misses only 0.122 MeV of 
the shell-model correlation energy) while the absolute energy of the ground state is rather 
worse underbound (by 0.407 MeV) in 6-IBM. 

The two-boson spectra obtained from mapping the four-hole shell-model results depend 
on the kind of pairs included in the basis but otherwise they are identical in the OAI 
and democratic mappings. Therefore, for each of the different IBM versions, there is a 
unique spectrum shown in Table IVl which is identical to that of the shell- model hamiltonian 
when diagonalized in the corresponding (possibly truncated) two-pair basis. While the 
OAI and democratic mappings yield the same energy spectrum for four holes, they lead to 
different boson-boson interaction matrix elements. Hence, the OAI and democratic results 
are different for the six- and eight-hole spectra shown in Tables I VI I and I VI II 

If the number of two-pair states equals the number of independent four-hole shell-model 
states, then the two-boson calculation reproduces the four-hole results exactly. According 
to Table IVl this happens, for example, for the three J'^ = 0^ states with T = which can be 
exactly described as combinations of |S'^;0), \D^]0), and \G^;0). Consequently, the three 
shell-model O"*" states are exactly reproduced in sdg-lBM. 

Because of the Pauli exclusion principle, no four-hole shell-model state exists with J = 18 
while this angular momentum is allowed in the coupling of two bosons with J = 9. This is an 
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TABLE VI: Energies (in MeV) of T = levels for six nucleons in the 1^9/2 orbit (^^Ag) calculated 
with the shell-model interaction SLGTO and compared with various versions of the IBM obtained 
with two methods of mapping, democratic (Dem) and OAI. Eq is the binding energy of the ground 
state. 
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example of the complication mentioned at the end of Sect. IIIII and the solution given there 
should be applied. In this case it implies that the matrix element (i?^; ISlHflB"^; 18) should 
be taken infinitely repulsive and it is only by adhering to this procedure that reasonable 
results can be obtained. 

Not much is known experimentally about ^^Ag except for the presence of two isomers, 
with tentative spin-parity assignments (presumably the lowest T = state) and 21"'" 
at about 5.78(3) MeV above the 7"*" [3]. The different shell-model interactions SLGTO, 
GF, and Nb90 all predict a 7"*" as the T = ground state, and a 21+ level at 6.464, 5.948, 
and 4.632 MeV, respectively. The 6-IBM reproduces the shell-model result for the binding 
energies of these isomers to about 100 keV for the 7"*" and less than that for the 21+ (see 
Table rvT] for the results of the SLGTO interaction). This result is valid for the different shell- 
model interactions: although the binding energies calculated with the three shell-model 
interactions vary by several MeV, in each case they are matched by the (appropriately 
mapped) 6-lBM to within about 100 keV, indicating that the B pair incorporates most of 
the correlations for the 7+ and 21"*" states. 

The 6-IBM results should be contrasted to those obtained with sd-IBM which fails com- 
pletely to reproduce the spectroscopy of ^^Ag. This is not surprising since it is known from 
the work of Elliott and Evans 29| that IBM-3 cannot give a satisfactory description of odd- 



odd nuclei which require the addition of isoscalar s and d bosons leading to the so-called 
IBM-4. While the latter is a realistic model when low-j shell-model orbits are involved {e.g., 



for S(i-shell nuclei 



i Q, 



311]), the present results seem to indicate that the mapping from a 



shell-model space with high-j orbits calls for the inclusion of an aligned isoscalar n-p pair 
with J = 2j. 

While 6-IBM and sb-lBM are largely equivalent for the odd-odd nucleus ^^Ag, this is 
not the case for the even-even nucleus ^^Pd. As can be seen from Table IVIIt the s boson 
provides crucial additional correlation energy which brings the boson result close to its shell- 
model equivalent. This nucleus was studied recently in a fusion-evaporation experiment j^. 
The excitation energies of the yrast levels, calculated with the SLGTO interaction, are in 
reasonable agreement with the observed values of 0.874, 1.786, and 2.535 MeV, respectively. 
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TABLE VII: Energies (in MeV) of T = levels for eight nucleons in the lgg/2 orbit (^^Pd) 
calculated with the shell-model interaction SLGTO and compared with various versions of the IBM 
obtained with two methods of mapping, democratic (Dem) and OAI. Eq is the binding energy of 
the ground state. 
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D. Electric quadrupole properties 

A further test of the aligned-n-p-pair hypothesis can be obtained from electric quadrupole 
(E2) transition properties. The E2 operator in the shell model is given by 

f^{E2) = ^ r^Y,^{ei, <f>i) + e^Yl ^l^^^iOi, <l>i), (46) 

where the sums are over neutrons and protons, and each sum is multiplied with the appro- 
priate effective charge. This operator can be written alternatively as a sum of an isoscalar 
operator, multiplied by {ci, + e^), and an isovector operator, multiplied by (e,^ — e^r). For 
the E2 transitions between T — levels of interest here, only the former part contributes. 
In second quantization, which is a convenient formalism for carrying out the mapping, the 
fermion E2 operator can be written as 

(47) 

where creates a neutron (p — u) ot a proton (p = tt) in the 1519/2 orbit, and djm — 
(— )-'+"*aj_TO. Furthermore, the factor in front comes from the radial integral over harmonic- 
oscillator wave functions (with length parameter Zho) involving the 15^9/2 orbit. 
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FIG. 2: Absolute values of the E2 reduced matrix elements for the transitions J ^ J — 2 between 
four-nucleon-hole states with T = in the 1(79/2 orbit (^^Cd), calculated with the SLGTO shell- 
model interaction and compared with the mapped 6-IBM. Matrix elements are expressed in units 
v/55/37r(et, + e^)/^^ (see text). 

The lowest-order bosonic image of the fermion E2 operator is defined by the diagonal 
(reduced) matrix element in the Q"*" state of the In-lp system which is given by 

{{lg,l,f-X\\T\m\{-^9s^/2?-X) = -^^IL X ^^(e. + e.). (48) 

The E2 operator of the 6-IBM is of the form 

f^{E2) = e^{b^xbY^\ (49) 

and is necessarily of scalar character in isospin. Since the mapping implies the equality 

((l(79/2)^9+||f^(E2)||(l(?9/2)^9+) = (6||T^(E2)||&), (50) 

and the since the boson matrix element on the right-hand side equals y/be^, we find the 
following expression of the boson effective charge Cb in terms of the shell-model neutron and 
proton effective charges: 




In the following, the factor y 55/37r(e,y+e,r)^ho is divided out of all matrix elements, fermionic 
as well bosonic. 

A first test is from E2 transitions between four-nucleon-hole states with T = 0. The shell- 
model results obtained with the SLGTO interaction, shown in Fig. |2l display a characteristic 
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FIG. 3: Absolute values of the E2 reduced matrix elements for the transitions J — )• J — 2 between 
six-nucleon-hole states with T = in the 1^9/2 orbit (^^Ag), calculated with the SLGTO shell- 
model interaction and compared with the mapped 6-IBM. Matrix elements are expressed in units 
v/55/37r(et, + e^)/^^ (see text). 
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FIG. 4: Same caption as Fig. [3] for J — )• J — 1 transitions. 

decrease in quadrupole strength for J ^ 8 which can be viewed as a remnant of a seniority- 
like classification. The figure also shows the results found in 6-IBM using the boson effective 
charge derived in Eq. (15T]) . with no adjustable parameter. Not surprisingly, given that the 
J = 8 state is poorly described in terms of aligned n-p pairs (and hence b bosons), the two 
transitions involving this state deviate strongly in 6-IBM from the corresponding shell-model 
result. Other transitions agree, however. 

A second test is provided by E2 transitions between six-nucleon-hole states with T = 0. 
They are shown in Figs. |3]and|l]for J — )■ J — 2 and J — J — 1, respectively. Assuming 
that an agreement between shell model and 6-IBM is obtained only if both the initial and 
final states are adequately represented by b bosons, we conclude that the 6-IBM is a good 
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FIG. 5: Same caption as Fig. [2] for eight-nucleon-hole states (^^Pd). 

approximation for two ranges of angular momenta, namely J = 6 to 13 and J = 17 to 21. 
This conclusion agrees, at least qualitatively, with the one drawn on the basis of energies (see 
Table rvTI|) . Since three b bosons cannot couple to total angular momentum J = 2, this state 
is absent from 6-IBM while it is present in the shell model (see Table IVlT|) . As a consequence, 
no 4 — 7- 2 or 3 —> 2 transitions occur in 6-IBM (see Figs. [3]and|l]). These transitions exist in 
the shell model but it is rather striking that the calculated matrix elements are very small 
indeed. 

Finally, in Fig. |5] are shown the E2 transitions between eight-nucleon-hole states with 
T = 0. A small depletion of the E2 strength calculated in the shell model is perceptible for 
J ^ 8 and is absent in 6-IBM. Apart from this deviation both calculations agree, indicating 
that the shell-model wave functions can be adequately represented in terms of a single b 
boson. We emphasize once more that, the number of states that can be written in terms 
of B pairs is small (of the order of 10), no effective charges are needed to arrive at the 
agreement in Fig. |5l 

Let us now formulate an overall evaluation of these results. While the 6-IBM gives in 
all cases an reasonable description of the ground-state binding energy, the addition of the 
standard s boson (with J = 0,T = 1) further improves the agreement. In fact, the energies 
obtained in sb-lBM (both with the OAI and democratic mappings) agree well with those 
of the shell- model levels except for (i) the J = 8 level in the four-hole system, (ii) low- J 
levels of the six-hole system, and (iii) levels with J = 14 to 16 in the six-hole system. The 
first discrepancy is obviously related to the small overlap of the J = 8 shell-model state 
with J = 8), noted in Table HVl and explained in Sect. IIVBI The second difference is 
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also understandable since a correct description of the low- J states in the odd-odd nucleus 
requires the consideration of low- J pairs with T = which have been omitted from the 
present mapping. The third deviation could be related to an unfavorable coupling of three 
B pairs to the angular momenta J = 14, 15, and 16, akin to the coupling of two B pairs to 
J = 8. The results as regards E2 transitions are consistent with what is concluded from the 
analysis of spectra. 



V. CONCLUSIONS 



We have shown in this paper that part of the low-energy spectroscopy of N = Z nuclei 
which have their valence nucleons confined to a single high-j orbit, can be represented in 
terms of an aligned isoscalar n-p B pair with J = 2j and is further improved by the inclusion 
of an isovector 5* pair coupled to J = 0. This was proven explicitly for a four-hole system 
and indirectly, through a mapping onto a corresponding boson system, for six and eight 
holes. Some deficiencies were found in this approach. A first concerns states of the four-hole 
system with angular momentum J ^ 2j which turn out to be poorly approximated with 
just S and B pairs. A second deficiency is more generic and pertains to the low- J states in 
odd-odd N = Z nuclei, the description of which calls for the inclusion of isoscalar n-p pairs 
with low angular momentum. Nevertheless, it should be noted that the two isomers that 
have been observed so far in ^"^Ag, (T"*") and (21+), are adequately described in terms of b 
bosons. 

These results were obtained for the lg9/2 orbit and for three different choices of two- 
body interaction. To what extent are they valid in general and can they be considered as 
representative of a system of neutrons and protons confined to a high-j orbit? In essence, 
two ingredients, geometry and dynamics, determine the outcome of the present pair analysis. 
The geometry is defined by the CFPs and, provided j is not too small, is expected to evolve 
only slowly with j. (It would in fact be an exercise of some interest to perform the pair 
analysis in the limit of large j.) The dynamics is determined by the two-body interaction 
which in our study was varied significantly but within reasonable bounds. The matrix 
elements shown in Table [I] are typical of what is obtained for a residual interaction with a 



short-range, attractive character 
to orbits other than lgg/2- 



32l | and we may thus expect similar results when we move 
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This work calls for further studies. The pair analysis of the shell-model wave functions 
should be extended to higher particle numbers which can be achieved through an isospin- 
invariant formulation of the nucleon-pair shell model. Consequently, the present results 
require further confirmation at higher particle number but one is tempted to conclude at 
this point that a realistic model can be formulated in terms of s and b bosons. Due to its 
simplicity, such a model could be of use to elucidate the main structural features of N Z 
nuclei in this mass region. These topics are currently under study. 
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